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1 Introduction 

As this volume indicates, the technology of atom chips is currently enjoying great success for a 
large variety of experiments on degenerate quantum gases. Because of their geometry and their 
ability to create highly confining potentials, they are particularly well adapted to realizing one 
dimensional (ID) situations [1-9]. This characteristic has contributed to a revival of interest in 
the study of ID Bose gases with repulsive interactions, a system which provides a vivid example 
of an exactly solvable quantum many-body system [10-12]. The quantum many-body eigenstates 
[10, 11] and thermodynamics [12] can be calculated without resorting to approximations. In 
addition, the ID Bose gas shows a remarkably rich variety of physical regimes (see Fig. [D 
that are very different both from those found in 2D and in 3D. One dramatic example of the 
difference is the tendency for a ID Bose gas to become more strongly interacting as its density 
decreases [10]. Finally, and in a more practical vein, a good understanding of its behavior is 
relevant for guided-wave atom lasers [13] and trapped-atom interferometry [14]. Because of the 
effects of interactions, the analogy to the manipulation of light in single mode fibers needs to be 
examined carefully. 

An atom chip is not the only means of producing a ID Bose gas. Optical trapping has been 
used to generate similarly elongated trap geometries. In particular, a 2D optical lattice can be 
used to generate a 2D array of ID traps [15-19]. Because of the massively parallel nature of this 
system, it is possible to work with only a few atoms per tube, and still get a sizeable signal per 
experimental cycle. Thus, the strongly interacting regime alluded to above can be reached. This 
regime has yet to be reached with an atom chip. But as we will show here, a key feature of atom 
chips is that they produce individual samples in which one does not intrinsically average over 
many realizations. Fluctuation phenomena are therefore readily accessible, an aspect which we 
will treat later in this chapter. 

In the following we first give an introduction to the various regimes of the homogeneous 
ID Bose gas, with particular emphasis on the behavior of the density profiles and the density 
fluctuations in the context of approximate models. Then we will discuss the exact solution and 
how it differs from the approximations. Next, we discuss some of the important issues involved 
in realizing ID gases in a 3D trap. Finally, we describe a series of experiments performed in 
Orsay and Amsterdam using atom chips to explore and illustrate features of the ID Bose gas. 

2 Regimes of one-dimensional gases 

First, we review some theoretical results concerning the one-dimensional Bose gas with repulsive 
interactions. Most of these results are derived in Refs. [10, 12,20-24]. Here we will concentrate 
on intuitive arguments, and the reader is referred to the above references for more careful demon- 



2 



strations. The system is described by the Hamiltonian 



H = / dzip + —ip + - I dzip + ?p + i;ip, (1) 



2m J d 2 J 2 

where ip is the field operator in second quantization, and g is the coupling constant characterizing 
the interactions between particles. From this coupling constant, one can deduce an intrinsic 
length scale related to the interactions, 

h 2 

l a = — , (2) 
mg 

as well as an energy scale: 

mg 2 h 2 

9 = 2fi? = 2m7f " (3) 

In thermal equilibrium, the gas is described by the temperature T and the linear atomic den- 
sity n. Rescaling these two quantities by the intrinsic scales introduced above, and setting Boltz- 
mann's constant equal to unity (i.e., measuring temperature in units of energy) we find that the 
properties of the gas are functions of the dimensionless quantities 

T 

and 

= mg_ = _\_ 

h 2 n nl g ' 

the latter being the famous Lieb-Liniger parameter [10]. 

It is useful to also introduce two other relevant scales, namely the thermal de Broglie wave- 
length, 

/ 27T [Ati 

^dB = fl\ —=; = L\ —, (6) 



and the quantum degeneracy temperature 



mT 



h 2 n 2 _ E g 
2m 7 2 



(7) 



In the above (t, 7) parametrization, quantum degeneracy (T w T d , or equivalently n\ dB ps 1) is 
reached around 

t ~ \- (8) 

T 

The thermal equilibrium for the hamiltonian of Eq. ([T]) has been extensively studied theoreti- 
cally [12,22]. Without going into great detail however, we can present some important features of 
this system. Several regimes may be identified in the parameter space (7, t), as sketched in Fig.[TJ 
We begin by noting that the region 7 > 1, t < 1 (dark grey area) defines a strongly interacting 
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Figure 1 : Physical regimes of a ID Bose gas with repulsive contact interactions in the parameter 
space (7, t), adapted from [22]. The dashed diagonal line separates the degenerate and nonde- 
generate gases. The strongly interacting regime is shown in dark grey. The weakly interacting 
regime is divided into the nearly ideal gas regime (also called decoherent regime) shown in white 
and the quasicondensate regime shown in light grey. Note that the nearly ideal gas can be de- 
generate. The quasicondensate regime is divided into the thermal and quantum regimes. The 
lines represent smooth ( and often wide ) crossovers rather than phase transitions. The crossovers 
are given in Eqs. rfiOl) . rfTTT) . ((§]), rf25l) and rt471) . The dashed area shows the parameter space 
investigated in the experiments presented in this chapter. 



regime that occurs at low density and low temperature, often referred to as the Tonks-Girardeau 
gas [20,25,26]. 

In the weakly interacting regime, 7 < 1, several sub-regimes are identified. These are the 
regimes which to date have been accessible in atom chip experiments, and we shall elaborate 
further on their nature in the discussion below. The two main regimes are the nearly ideal gas 
regime (white area) and the quasi-condensate regime (light grey area). Each one permits an 
approximate description that we present later in this section and which allows the identification 
of sub-regimes. For the moment we simply wish to emphasize that no phase transition occurs 
in the ID Bose gas and that all the boundaries represent smooth (and often broad) crossovers in 
behavior. 

2.1 Strongly versus weakly interacting regimes 

We first comment on the distinction between strong and weak interactions. Following the ap- 
proach of Ref. [20], we study the scattering wave function of two atoms interacting via the 
potential g8{z\ — z 2 ), where z\ and z 2 are the position of the two atoms. For this, we consider the 
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wave function ip in the center-of-mass frame, with reduced mass m/2 and subject to the potential 
gS(z). The effect of the potential is described by the continuity condition 

JV( 0+ ) - |,xo_) = ^(o) (9) 

where + (0_) denotes the limit when z goes to zero through positive (negative) values. Let us 
consider the scattering solution for an energy E = h 2 k 2 /m. Since we consider bosons, we look 
for even wave functions of the form cos(/c|z| + (p). The continuity conditions give <\> and thus 
the value ^(0). We find then that the energy E g given by Eq. © is the relevant energy scale and 
that for E <C E g , ip(0) is close to zero, while, for E 3> E g , tp(0) is close to one, as illustrated in 

Fig.m 

The above results hold for a gas of particles since the continuity relation © holds for the 
many-body wavefunction when two atoms are close to the same place. Thus, as long as the 
typical energy of the particles is much lower than E g , the many-body wavefunction vanishes 
when two particles are at the same position: the gas is then in the strongly interacting, or Tonks- 
Girardeau regime. The vanishing of the wave function when two particles are at the same place 
mimics the Pauli exclusion principle and the gas acquires some similarities with a gas of non 
interacting fermions. More precisely, in this strong interaction regime, the available wave func- 
tions of the many body problem are, up to a symmetrization factor, the wave functions of an 
ideal Fermi gas [26]. Since the wave function vanishes when two atoms are at the same place, 
the energy of the system is purely kinetic energy and the eigen energies are those of the Fermi 
system. Thus the ID strongly interacting Bose gas and the ideal ID Fermi gas share the same 
energy spectrum. This implies in particular that all thermodynamic quantities are identical for 
both systems. 

To identify the parameter space of the strongly interacting regime, we suppose the gas to be 
strongly interacting and then require that the typical energy of the atoms be smaller than E g . To 
estimate the typical energy per atom, we use the Bose-Fermi mapping presented above. If the 
gas is degenerate, the temperature is smaller than the degeneracy temperature Td, Eq. ©, and Td 
corresponds to the "Fermi" energy of the atoms. The typical atom energy is therefore Td and it 
is of order E g if 

7 ~ 1. (10) 

The strongly interacting regime thus requires 7 3> 1. If the gas is non degenerate, the typical 
energy of the equivalent Fermi gas is T and interactions become strong when T = E g or 

f~l. (11) 

We then find that the gas is strongly interacting for t < 1. 

The condition (flOl) is often derived using the following alternative argument, valid at zero 
temperature. At zero temperature, there are two extremes for the possible solutions for the wave 
function ip(zi,Z2, •••)• As seen in Fig. |2l either the wave function vanishes when two atoms 
are at the same place, or the wave function is almost uniform, corresponding to the strongly 
and weakly interacting configurations respectively. In the weakly interacting configuration, the 
kinetic energy is negligible and the interaction energy per particle, of the order of gn, determines 
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Figure 2: Strong interaction versus weak interaction regime. We show the wave function in the 
center-of-mass frame of two atoms for (a) strong interactions, scattering energy E much smaller 
than Eg = mg 2 /2h 2 and (b) weak interactions, E much larger than E g . We also plot the wave 
function ip(z ly z 2 , z 3 , ...) for given positions ofz 2 , z 3 , ... in (c) the strongly interacting regime and 
(d) the weakly interacting regime. 

the total energy. In the strongly interacting configuration, on the other hand, the interaction 
energy vanishes while the typical kinetic energy per particle is h 2 n 2 /m. Comparing these two 
energies, we find that the strongly interacting configuration is favorable only for 7 > 1. 

2.2 Nearly ideal gas regime 

At sufficiently high temperatures, interactions between atoms have little effect and the gas is 
well described by an ideal Bose gas. In Ref. [22], this regime was referred to as the "decoherent 
regime"; We will call it the (nearly) ideal Bose gas regime. A ID ideal Bose gas at thermal equi- 
librium is well described using the grand canonical ensemble, introducing the chemical potential 
fi. All properties of the gas are calculated using the Boltzmann law which states that, for a given 
one-particle state of momentum hk, the probability to find N atoms in this state is proportional 
to e~^ 2fc2 li 2m )-f J ) N / T • no t e that \i < in this description. In the following, we use a quantization 
box of size L (tending to infinity in the thermodynamic limit) and periodic boundary conditions 
so that the available states are the momentum states with momentum k = 2irj/L where j is an 
integer. 

Let us first consider the linear gas density. From the Boltzmann law, we find that the mean 
population is the Bose distribution 

V 1 * 1 / = e (h 2 k 2 /(2m)-fj,)/T _ I ' 

The atom number, and thus the linear density, is obtained by summing the population over the 
states and one finds 

n=^-g 1/2 ^ T ), (13) 

*dB 

where g x / 2 (x) is one of the Bose functions 

i=i 
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also known as the polylogarithmic functions [27, 28]. Unlike in 3D systems, where the excited- 
state density is given by p e = 5'3/2(e M//T )/A^ B in this approach [27], no saturation of the excited 
states occurs (the function gi/ 2 diverges as a/— vrT/V as yU — > from below, whereas #3/2(1) = 
2.612 is finite): in the thermodynamic limit no Bose-Einstein condensation is expected and the 
gas is well described by a thermal gas at any density. 

Two asymptotic regimes may be identified: the non degenerate regime for which —fj, 3> T 
and h 2 n 2 /m <C T and the degenerate regime for which — ji <C T and h 2 n 2 jm 3> T. In the non 
degenerate regime, the linear density is well approximated by the Maxwell-Boltzmann formula 

n = -!-e^ T , (15) 

In this regime n\dB is much smaller than unity. In the degenerate regime, the states of energy 
much smaller than T are highly occupied and the linear density is given by 



-W5 

This density is much larger than 1/XdB, i-e. n\ dB 1. 

As we will discuss in the experimental section, fluctuations are also very important for char- 
acterizing the gas. It is thus instructive to consider the correlation functions. The normalized 
one body correlation function is g^\z) = (^ + (0)^(z))/ra, where tp is the field operator in 
the second quantization picture. Using the expansion of the field operator in the plane wave 
basis ip(z) = Ylk a k e ~ lkz /VL where is the annihilation operator for the mode k, we find 
g^(z) = J2k( n k) e ~ lkz /(Ln). Here n k = a^a k is the atom number operator for the mode k. 
Simple analytical expressions are found in the nondegenerate and highly degenerate limits. In 
the non degenerate limit (— fi 3> T or, equivalently n<l/ A^s), we find 



gW(z) ~ e (17) 

As the gas becomes more degenerate, the correlation length increases and, in the degenerate 
regime (— ji <C T or, equivalently n 3> 1/A^b), we find 

„ 2-kz 
ml z 



g {1 \z) ~ = e . (18) 

In this regime the correlation length, about n\ 2 dB , is much larger than the de Broglie wavelength 
(and the mean interparticle distance 1/n) since XdB ^> 1/n. 

Next we consider the normalized density-density or two body correlation function 

g&>(z) = (tp + (z)ip + (0)^(0)^(z))/n 2 . (19) 

This function is proportional to the probability of finding an atom at position z and at position 
z = 0. It is given by 

n 2 g^(z)= £ ««iM^-*«/L a . (20) 

kik2k:j k± 
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Using Bose commutation relations and the fact that, since atoms do not interact, different mo- 
mentum state populations are uncorrelated, the sum simplifies to: 



n 2 g (2 \z) = (n kl )(n k2 )(l + e^~ k ^)/L 2 + J>+a+a fc a fc > / ' L 2 • (21) 

fci^fc2 k 

In the last term, the commutation relations give: (a^a^a k a k ) = (n 2 ,) — (rik), and in thermal 
equilibrium one has: 

(nl) = (n k ) + 2(n k ) 2 . (22) 

Therefore we find: 

g®(z) = l + \gV>(z)\ 2 , (23) 



a result which one can also obtain directly from Wick's theorem [29]. Equation (1231) means that 
the probability of finding atoms within less than a correlation length in a thermal Bose gas is 
twice that of finding two atoms far apart. This phenomenon is often referred to as "bunching" 
and has been observed in cold atoms in several experiments [30-32]. Bunching is closely related 
to density fluctuations. As one can see from Eq. (1221) . in a thermal gas, fluctuations in the occu- 
pation of a single quantum state, 5n k 2 = (nl) — (n k ) 2 , show a "shot noise" term, (n k ) and an 
"excess noise" term, (n k ) 2 . The density fluctuation experiment described later in this chapter has 
demonstrated this behavior. 

Validity of the ideal gas treatment. The two body correlation function has been used to 
characterize the crossover between the ideal gas and quasi-condensate regimes [22]. When in- 
teractions become important, they impose an energy cost on density fluctuations and the latter 
tend to smooth out. This amounts to a reduction in the value of g( 2 \0). In the quasi-condensate 
regime which we discuss in the next section, the bunching effect is absent and g^(0) is close 
to unity. The ideal Bose gas description fails when the typical interaction energy per particle gn 
is not negligible compared to — /i. Using Eq. ([TBI one finds that the ideal Bose gas description 
fails when the temperature is no longer much smaller than the crossover temperature, which we 
define as 

T co ~ T dy fy. (24) 
Using the reduced dimensionless temperature t = T/E g , this can be written as 

This line separates the nearly ideal gas regime from the quasi-condensate regime in Fig.[Q Note 
that, in terms of chemical potential, the domain of validity of the ideal gas model is —fi ^> a co 
where we define the crossover chemical potential as 

T 

Hco = (26) 

In making this estimate, we have assumed that the gas is degenerate at the crossover. From 
Eq. (1241) . one can see that if one is in the weakly interacting regime (7 1) this assumption 
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is indeed true. The experiments described below confirm that one can observe the effects of 
degeneracy before the onset of the reduction of density fluctuations. 

A precursor of the reduction of density fluctations is shown by a perturbative calculation valid 
in the nearly ideal gas regime which gives, to lowest order in g [22], 

gM(0) ~ 2 - 4(T co /T) 2 . (27) 

To accurately treat the crossover regime however, it is necessary to make use of the exact solution 
to the ID Bose gas model. The exact solution in the crossover regime is discussed in Sec. 12.41 

The correlation lengths of the gas are important parameters of the gas that will be used in the 
following to estimate the validity criteria of the local density approximation. In the degenerate 
regime, the correlation length is l c ~ n\ 2 dB (see Eq. ([781)). Using Eq. (|24l) . we find that, close to 
the crossover, the correlation length of the gas is close to the healing length 

t = -4=- (28) 



2.3 Quasi-condensate regime 

On the other side of the crossover, i.e. for T <C T co , the bunching effect is entirely suppressed 
and the function is close to unity for any z. This regime is the quasi-condensate regimeQ. 
In this section, we present a description of the gas, valid in the quasi-condensate regime. This 
description permits a simple estimate of the density fluctuations. We thus verify a posteriori that 
the quasi-condensate regime is obtained for T <$C T co . We also give a simple calculation of phase 
fluctuations in the quasi-condensate regime. 

In the quasi-condensate regime density fluctuations are strongly reduced compared to their 
value in an ideal Bose gas where the bunching effect is responsible for density fluctuations of the 
order of n 2 . In other words: 

5n 2 < n 2 (29) 

In this regime, a suitable description is realized by writing the field operator as ip = e %e \/n + 5n 
where the real number n is the mean density and the operator 5n and the phase operator 9 are 
conjugate: [6n(z), 0(z')] = iS(z — z'). Note that the definition of a local phase operator is subtle 
and the condition Eq. (|29l) is not well defined since, because of shot noise, 5n 2 is expected to 
diverge in a small volume. A rigorous and simple approach consists in discretizing the space so 
that in each cell a large number of atoms is present while the discretisation step is much smaller 
than the correlation length of density and phase fluctuations [33]. 

Following this prescription, one first minimizes the grand canonical Hamiltonian H — fxN 
with respect to n to obtain the equation of state 

\i = gn. (30) 

'it is also called coherent regime since the g( 2 > function is close to unity, as in a coherent state. On the other 
hand, the first order correlation function still decays and so the gas is not strictly coherent in this sense. Within this 
terminology, the ideal Bose gas regime is called the decoherent regime [22]. 
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To second order in Sn, this is the correct expression of the chemical potential. This equality 
ensures that the Hamiltonian has no linear terms in Sn and V6>. Linearizing the Heisenberg 
equations of motion in Sn and VO, we obtain [33] 

HdSn/dt = 2^/h~(-^A)9^n~ 

These equations are the so-called hydrodynamic equations. They are derived from a Hamiltonian 
quadratic in Sn and V6>, that can be diagonalized using the Bogoliubov procedure [33]. It is not 
the purpose of this chapter to detail this calculation and to give exact results within this theory. 
We will simply give arguments that enable an estimate of the density fluctuations and of their 
correlation length. This estimate will then be used to check that Sn 2 <C n 2 , as assumed in 
Eq. ( |29l ). We will show that this condition is the same as the condition T <C T co where T co given 
in Eq. (1241) . After that, we will give similar arguments to estimate the phase fluctuations. Since 
in the following we will study the gas properties versus the chemical potential, it is instructive to 
rewrite the condition T <C T co in terms of chemical potential. Using Eq. (l30l) . we find that the 
quasi-condensate regime is valid as long as u ^> fi co where n co is given by Eq. (l26l) . 



2.3.1 Density fluctuations 

To estimate the density fluctuations introduced by the excitations, it is convenient to divide the 
excitations in two groups: the excitations of low wave vector for which the phase representation 
is most appropriate and the excitations of high wave vector for which a particle point of view is 
most convenient. 

In the following, we use the expansions on sinusoidal modes 9 = J2k>o V2(9 ck cos(kz) + 
9 sk sm(kz)) and Sn = J2 k>0 \/2(Sn ck cos(kz) + 5n sk sm(kz)). Here 5rij k and 9 jk are conjugate 
variables ({5nj k ,9j> k /} = {i/ L)8jj>5 k y) where j stands for c or s. For modes of small wave 
vector k, the excitations are phonons, or density waves, for which the relative density modulation 
amplitude 5rij k / n is much smaller than the phase modulation amplitude 9j k . In this case, the local 
velocity of the gas is given by KV9 jm and the kinetic energy term is simply Lnh 2 k 2 9 2 k / (2m). 
The Hamiltonian for this mode then reduces to 

H jk = L (g5n 2 k /2 + nh 2 k 2 9 2 k /(2m)) . (32) 

This hamiltonian could also have been derived from the equations of motion given in Eq. (T3TI) . 
provided that the quantum pressure term h 2 / (2m)A5n/n is neglected: indeed, for a given wave 
vector k, the laplacians in Eq. (|3T| ) give a factor k 2 and Eqs. (|3Tb are simply the equations 
of motion derived from the Hamiltonian Eq. (|32i For temperatures much larger than ng, the 
thermal population of these phonon modes is large and classical statistics apply. Thus, the mean 
energy per quadratic degree of freedom is T/2 and we obtain 

(5n 2 k )=T/(Lg). (33) 

and 

(9 2 k ) = mT/(Lnk 2 H 2 ). (34) 
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We can now check the validity of the assumption 5rik/n <C 9 k : it is valid as long as k <C 
yjmgn / h. Since k values are spaced by 2n/L, there are about L^Jmgnj (irk) modes that satisfy 
this condition. Since the contribution of each of these modes to the relative density fluctuations is 
given in Eq. (l33l) . we find that the contribution of these low momentum excitations to the relative 
density fluctuations is of the order of 

(^ 2 )phonons _ T _ _T_ _ T/T ^ (35) 



nh-J 'gn/m Tdy/j 



For wave vectors much larger than yjmgn / h, the phase-density representation is not the 
most appropriate. An excitation of wave vector k ^> ^Jmgn / h corresponds to the presence of 
an atom of momentum k, whose wave function is e lkz / 'y/L and whose energy is h 2 k 2 /(2m). 
The anihilation operator for this mode is a k as introduced in section [2T2] For temperatures much 
larger than h 2 k 2 /m, the thermal population of this mode is large and classical field theory, in 
which a k is treated as a c-number, is adequate. We then find that a k has a Gaussian distribution 
which satifies (|ctfc| 2 ) = 2mT / [h 2 k 2 ). The density fluctuations caused by the presence of such 
high momentum atoms result mainly from the interference between the atomic field a k e %kz j\fTj 
and the atomic field of long wavelength spatial variations, whose amplitude is close to ^/n. The 
density fluctuations are thus 5n = \/n(ake lkz + a* k e~ lkz ) / \TL. We then find that the contribution 
of the mode of wave vector k to density fluctuations is Sn 2 , = AnmT '/ '(Lh 2 k 2 ). Summing the 
contributions of the modes for all k > ^Jmgn/ h, we obtain an estimate of the density fluctuations 
(<fo 2 ) atoms caused by high momentum excitations: 

(Sn 2 ) atoms T T 

2 - ™ ^ - ~- ( 36 ) 



n 2 T dv /7 T t 



CO 



One also sees from the above argument that the density fluctuations fall off as 1 / k 2 above k = 
y/mgn/ h. The inverse of this scale gives the length scale of density fluctuations and we find that 
this correlation length is the healing length £ defined in Eq. (|28l) . 

From Eq. (1351) and Eq. (l36l) . we find that Sn 2 /n 2 ~ T/T co . Thus, the quasi-condensate 
treatment is valid as long as T T co , In conclusion, we have shown that T co gives the limit 
of both the ideal gas regime, valid as long as T 3> T co , and the limit of the quasi-condensate 
regime, valid for T T co . Equivalently, in terms of chemical potential, as long as the chemical 
potential is positive and much larger than /i co of Eq. ( |26l ), the gas is in the quasi-condensate 
regime whereas for negative chemical potential of absolute value much larger than /i co the gas 
is in the ideal gas regime. This is illustrated in Fig. [31 The two regimes differ by the fact that 
the g^'{z) function is modified: it is close to one for any z in the quasi-condensate regime while 
</ 2 )(0) = 2 in the ideal gas regime. 



2.3.2 Phase fluctuations 

In the quasi-condensate regime, although the gas is coherent with respect to the g^ function, it 
is not coherent with respect to the g^ function. This is why the gas is called a gwasz-condensate. 
The phase fluctuations have been measured experimentally in various experiments where the 
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quasi-condensate presented a one-dimensional character [34-38]. The description of the quasi- 
condensate given above permits a simple calculation of those phase fluctuations as we now show. 
Phase fluctuations are given by 

((6(z) - 6(0)) 2 } = ^2(^)(cos(fe) - l) 2 + 5><0 s 2 fe > sin 2 (kz). (37) 

fc>0 fc>0 

Using Eq. (l34l) and (cos(fcz) — l) 2 + sin 2 (/cz) = 2(1 — cos(kz)) this gives 

((9(z) - 9(0))*) = A(mT/(Lnh 2 )) £ L^M . (38 ) 

fc>0 

Transforming into L/ (27r) J 00 dfc and using J °°(l — cos(kz))/k 2 dk = irz/2, we obtain 

((0(z)-0(0)) 2 } = ^ = ^. 09) 
nh 2 n\ 2 dB 

Since density fluctuations are very small, the function is about g^ l \z) = n(e^ e ^~ 9 ^). 
Since the Hamiltonian is quadratic, we can use the Wick theorem to compute (e^ 6 *^ -0 ^), 
which gives ( e *(*(*M(o))) = e -<(^M(o)) 2 >/2_ We find 

gM(z) ~ e - mTz/{2nh2) . (40) 

Comparing this to Eq. (fT8l . we observe that the behavior of g^ 1 ' is close to that in the ideal gas 
regime. The factor of 2 difference in the correlation length formulae is because for the ideal 
gas regime, both density and phase fluctuations contribute to g^ whereas only phase fluctua- 
tions remain in the quasi-condensate regime. The crossover from the ideal gas regime to the 
quasi-condensate regime, at a temperature T co , Eq. (|25T) . corresponds to the situation where the 
correlation length of phase fluctuations, given by Eq. (1401) . equals the correlation length of den- 
sity fluctuations given by Eq. (|28l ). 

In both this section and the previous one, we assumed that the temperature is high enough that 
the population of the relevant modes (whose wavelengths are of the order of £) is much greater 
than unity. This is no longer the case when T reaches values of the order or smaller than gn. 
For lower temperatures, quantum fluctuations are expected to be dominant. This is the so-called 
quantum quasi-condensate and the boundary between the thermal quasi-condensate regime and 
the quantum quasi-condensate regime is at T w ng, corresponding to 

t»- (41) 

7 

and is shown as a line in Fig. CD A recent experiment using an atom chip observed these quantum 
phase fluctuations [39]. In the experiments we describe here however, the temperature is high 
enough that thermal fluctuations dominate. 
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2.4 Exact thermodynamics 



In sections l2T2l and l231 we have discussed models that apply independently in the asymptotic lim- 
its of the nearly ideal gas regime (T 3> T co or equivalently —fx 3> p co ) and the quasi-condensate 
regime (T <C T co or equivalently fi ^> p, co ) respectively. While the above classification gives 
very useful insight, it should be emphasized that the boundary between these two regimes is a 
smooth crossover, not a sharp transition and that neither of the two theories presented above ac- 
count for the physics in the vicinity of the crossover. Since in many cases we are interested in 
the precise behavior near the crossover from the ideal gas to the quasi-condensate regime, it is 
not sufficient to use the asymptotic results. 

As already mentioned in the introduction, the ID Bose gas with repulsive delta-function 
interactions is an example of an exactly solvable model [40,41], This allows us to quantita- 
tively compare predictions of the two approximate descriptions to the exact results, and verify 
the regions of validity of the approximations. Furthermore the exact results will turn out to be 
important for an accurate description of the experiments. 

Exactly solvable models typically occur in lower dimensions (ID quantum systems [40,41] 
and 2D classical systems [42]) and allow one to obtain exact solutions for the quantum many- 
body eigenstates through a method known as the "Bethe Ansatz" (due to Hans Bethe [43]), 
for any value of the interaction strength. For the repulsive delta-interacting ID Bose gas (with 
periodic boundary conditons), these solutions were first obtained by Eliot Lieb and Werner Lin- 
iger [10, 1 1]. Furthermore, the method based on the Bethe Ansatz can be extended to also obtain 
the thermodynamics exactly (for any temperature), via a method due to C. N. Yang and C. P. 
Yang [12]. 

For a concise and lucid description of the Yang-Yang method to obtain the exact thermody- 
namics of the ID Bose gas and the related equations, we refer the reader to the original litera- 
ture [12]. In brief, each exact quantum many -body eigenstate of the Lieb-Liniger hamiltonian 
Eq. CQ) is characterized by a set of distinct integer quantum numbers and a corresponding set 
of distinct quasi-momenta k, obtained through the Bethe Ansatz. For a large system, one can 
consider the distribution of these quasi-momenta p(k) and also of the "holes" ph{k), the latter 
corresponding to the "missing" values in the set of integers characterizing the individual quantum 
states. By considering the entropy for given distributions p(k) and ph{k), Yang and Yang showed 
that the condition of thermal equilibrium leads to a set of nonlinear integral equations that can 
be solved by iteration. Subsequently, from the resulting distributions thermodynamic quantities 
such as pressure and free energy can be obtained. Once these quantities have been found, further 
thermodynamic quantities can be calculated using the standard thermodynamic relations. 

Although numerical solutions to the Yang- Yang equations were already obtained at an early 
stage by C. P. Yang [44], important further insight into the Yang- Yang thermodynamics was 
gained much more recently by Kheruntsyan, Gangardt, Drummond and Shlyapnikov [22,23]. 
They calculated both density and the normalized local density-density correlation function g^> (0), 
and compared to approximate results in the various regimes discussed above. The former, 
n(p, T), is obtained as part of the equation of state. The latter is obtained from the derivative of 
free energy with respect to the coupling constant g, using the Hellmann-Feynman theorem. 

As an important example, a comparison to the approximate results of the previous sections 
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is shown in Fig. [3l for a fixed scaled temperature of t = 1000. This value is in the relevant 
range for the experiments to be described below. Such curves as a function of chemical potential 
fx are particularly useful to describe the behavior in a trap, since in this case one has a well- 
defined global temperature, while the density varies (within the local density approximation) 
with the local chemical potential fx(z) according to fx(z) — fx — V(z), where V(z) is the trapping 
potential. This will be discussed in more detail in the following section. 

Figure [3t a) shows that the exact density (n oc I/7) indeed approaches the ideal-gas behavior 
as [i I [i co becomes sufficiently negative, while for large positive fx/ fx co it approaches the quasi- 
condensate result. There is a large range in density (more than a factor 4) over which neither 
asymptotic description gives correct predictions. In the same vein, the local density-density 
correlation function g^(0) (Fig. Ob)) smoothly crosses over from 2, the value for an ideal gas, 
to about 1, as expected for a quasi-condensate. This smoothness is characteristic of crossover 
behavior, and is drastically different from the step-like behavior typical for a 3D gas. 

Looking more closely at Fig. [31 one sees that the ideal gas description begins to fail for a gas 
that is only moderately degenerate: already at fx/T = —0.5 (fx/ fx co = — 5 for the considered t 
parameter), a chemical potential for which nXas ~ 10 and the population in the k = mode 
according to Eq. (fT2l) is ps 1.5, the ideal Bose gas prediction is off by about 10%. This is because 
the interaction-induced crossover is sufficiently wide that for the used value of t (1000), the 
chemical potential at degeneracy (fx/T ~ —1) is not very far removed. The narrowness of the 
degenerate ideal gas regime is also seen in Fig.Q] To achieve well separated regimes, one would 
need to work at much higher t and much smaller 7. For t = 1000, the effect of degeneracy 
is nevertheless visible before the quasi-condensate crossover. This is shown by comparing the 
density with both the true ideal gas model and the Maxwell-Boltzmann model in Fig. 3J at 
fx/ fx co ~ —5 the ideal gas model gives a prediction for the density accurate within 10% as 
mentioned above (and the nearly-ideal-gas description can thus be expected be applicable) while 
the Maxwell-Boltzmann prediction is off by a factor of about 2. 

Concerning the local pair correlation function ^^(O), it deviates from the ideal-gas value of 
2 for the entire range plotted in the figure. The experiments presented in this chapter however 
(see Sec. |4]), are not precise enough to detect this deviation. Finally, the value of g^ 2 \0) can 
take values below unity in the quasi-condensate regime. We will briefly return to this point in 
Sec. 14331 

Despite its power, the Yang- Yang theory does not permit calculation of any non-thermodynamic 
quantities. For example, only the local value of the density correlation function g^(0) has be 
obtained from thermodynamics, while the full behavior of g^(z) has been obtained from the 
exact solution only at zero temperature [46]. At finite temperature, the behavior of g^ 2 \z) has 
been obtained only by perturbative calculations valid in each asymptotic regime [47], but they do 
not describe the crossover itself. An alternative approach uses the fact that the crossover appears 
in a highly degenerate gas. In this case, the modes are highly populated and a classical field 
approach is possible [48,49]. 
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Figure 3: Normalised density (I/7) and local pair correlation g^ 2 '(0) as a function of chemical 
potential scaled to the crossover chemical potential /x co given in Eq. ft26\) for fixed temperature 
corresponding to t = 1000. Numerical results from the Yang-Yang equations (solid lines, cour- 
tesy K. Kheruntsyan ) are compared to the ideal Bose gas model ( dotted line in (a), Eq. 4731) ) and 
the quasi-condensate model ( dashed line in (a), Eq. 001 ) ). The vertical line in (a) indicates the 
degeneracy chemical potential —\i = T. The classical Maxwell-Boltzmann prediction Eq. 4771 ) 
is shown as dashed-dotted line. In (b) the asymptotic values of g^(0) are indicated for both 
the ideal-gas regime (g^ 2 '(0) = 2 for p <C — p co , dotted line), and the quasi-condensate regime 
(g^(0) = lfor p ^> p co , dashed line). Adapted from Refs. [38,45]. 
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3 ID gases in the real world 



In real experimental situations, the atomic gas is neither homogeneous nor purely one-dimensional 
As usual, in our experiments the trapping is to a good approximation harmonic. The trap has 
cylindrical symmetry and is characterized by a tight radial trapping frequency u± and a much 
lower axial trapping frequency u. Here, we briefly summarize the main issues related to realiz- 
ing a ID system in this trapping geometry. We first discuss the link between transverse effects 
related to cu± and we present a model based on the Yang- Yang thermodynamics, valid at low 
enough linear densities, that takes into account these transverse degrees of freedom. We then 
discuss the effect of the longitudinal trapping potential. We finish by discussing the link with the 
3D physics, in particular with regard to the usual Bose-Einstein condensation in 3D. 

3.1 Transverse trapping and nearly ID Bose gases 

Strictly speaking the conditions to be ID in a transversely trapped gas are that both temperature 
and chemical potential are much smaller than the radial vibration quantum, T, [i <C fauj±. If this 
is the case, the gas is frozen in the transverse direction both thermally and in terms of chemical 
potential, and the (many-body) wave functions can be factorized into the product of a transverse 
part (the gaussian ground-state wavefunction of the radial trap) and an axial part. The system is 
then kinematically one-dimensional. Studying the scattering properties, Ref. [20] has shown that 
the interactions can be modeled by an effective ID coupling constant g and, as long as the 3D 
scattering length a is much smaller than the typical size of the transverse oscillator wavefunction, 
l± = \/K/nuJ±, 

g = 2aftw ± . (42) 

In most experiments on atom chips, neither of the above conditions on temperature and chem- 
ical potential are well-fullfilled, and is it necessary to also take into account the transverse degrees 
of freedom. 

It is useful to consider the linear density ni, obtained from the actual 3D density p(x, y, z) 
through integration 

ni(z) = J J dxdyp(x,y,z). (43) 

When the gas is strictly ID, one can identify rti with the ID density n. We will present our main 
experimental results in terms of this linear density, because it turns out that ni is often the key 
parameter, in particular when considering the crossover to 3D at low temperatures. 

This is in particular true for the quasi-condensate regime, at temperatures T <C fko±. In 
this regime, the chemical potential is close to its value at zero temperature, which is given by 
the solution of the radial Gross-Pitaevskii equation [50,51]. It was found from comparison to 
numerical integration of the radial Gross-Pitaevskii equation [50,5 1 ] that, in the quasi-condensate 



16 



regime, to good approximation the chemical potential can be expressed as^ 



fj, = hw± (VI + 4ra/a - l) . (44) 

For linear density rii <C l/4a, we find that /x ~ 2fku±an. We recover here the chemical potential 
gn of the ID case. At larger linear density, the chemical potential is reduced compared to the 
formula 2hcu±ani. This reflects the fact that, for large densities, the transverse cloud size is 
increased with respect to the transverse vibrational ground state. As another example of how n\ 
is the relevant quantity for low enough temperatures, we note that the expression Eq. (|4Q|) for the 
phase coherence length remains correct also on the 3D side of the crossover, if we replace the 
ID density n by the linear density n\ [52]. 

3.2 Applying ID thermodynamics to a 3D trapped gas 

Another case that one can consider is when the interaction energy is in the ID regime, ng <C fku±, 
while temperature is in the 1D-3D crossover, T ~ Tko±. A model for this regime was introduced 
in Ref. [8], and we describe it here. The key step is to separately consider the radial states. Under 
the above conditions only the radial ground state is significantly affected by the interactions, 
while the radially excited states can still be treated as an ideal gas. Thus, for the radial ground 
state, the solution uyy{h,T) to the Yang-Yang equations must be used. Each radially excited 
state with radial quantum number j > 1 is now considered as an independent ideal ID gas, in 
thermal equilibrium with the rest of the cloud. Each of the radially excited states is thus taken to 
have a density (cf. Eq. (fT3l) ) 

n e (fij,T) = —g 1/2 {exp(jJLj/T)), (45) 

where an effective chemical potential fij has been introduced that takes into account the radial 
excitation energy, 

Hj = fj, - jhw±. (46) 

Taking into account the degeneracy factor j + 1 of the radially excited states, the total linear 
density in this model thus becomes 

oo 

ni(n,T) = n YY (iJ,,T) + ^(j + l)n e (/jj,T). (47) 

3=1 

As long as /i < fku±, we have jij < which is necessary to avoid divergence of gyi in Eq. (|45T) . 
In fact, from the previous discussion in Sec. I2.4L for our parameters (t ~ 1000), we can expect 
an ideal-gas treatment of the radially excited density to begin to break down for fij/T > —0.5 
since this is where interactions will become important. In practical cases where T w hu ± , the 
model should thus be accurate as long as [i < 0.5fku±, while for ji > 0.5hw± the model will start 
to become inaccurate. 

2 An additional factor —1 has been introduced in brackets in Eq. d44l i compared to Ref. [51]. This subtracts the 
radial zero-point energy hu>±, so that ii = corresponds to the energy of the lowest energy (fc = 0) state, as in the 
treatment in Sec. [2] 
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3.3 Longitudinal trapping 

Experimentally, cold gases are axially confined in a confining potential V(z) and the cloud is not 
infinite and homogeneous as assumed in the previous section. However, as seen below, for weak 
enough axial confinement, the results for homogeneous gases can be applied using a local density 
approximation. In the first following sub-section, we present the local density approximation and 
discuss its predictions. We then evaluate the condition of validity of this approximation. 

3.3.1 Local density approximation 

If the confinement is weak enough that the correlation length of the gas is, at each position, much 
smaller than the length of the mean density variations, then the gas may be divided into small 
slices in which the thermodynamics of uniform systems applies. A slice located at position z is 
in equilibrium with the rest of the gas. It is thus described by the grand canonical ensemble at 
temperature T and at a chemical potential /i . The energy of the gas contained in this slice is 
shifted by the quantity V(z). It is equivalent to assuming that the chemical potential is fi — V(z), 
while the energy of the gas is unshifted. Thus, the local properties of the gas are that of a 
homogeneous infinite gas at temperature T and local chemical potential /jl(z) = /x — V(z). This 
is the so-called local density approximation. 

Within the local density approximation, all the results presented in the previous section hold. 
Thus, performing local analysis, one can observe all the features of homogeneous ID gases: the 
presence of the ideal gas regime, which includes the degenerate regime, the crossover towards a 
quasi-condensate and the quasi-condensate regime. In particular, a quasi-condensate appears in 
the center of the trap, when the peak density exceeds the crossover density n co given by Eq. (|24|) . 

It is often interesting to investigate the behavior of the gas using the extensive variable N, 
where N is the total atom number. As long as the peak density is much smaller than n co , the 
density profile is well described using the equation of state n(fx, T) of an ideal Bose gas. Then 
the total atom number is easily computed and, for gases that are degenerate at the trap center, we 
obtain [45] 

N = T/(huj)\n(T/\fi \). (48) 

The atom number at the crossover towards a quasi-condensate is obtained when the peak density 
reaches n co . Inserting Eq. (fT6l) and Eq. (1241) into Eq. (|48l) we find that the atom number at the 
crossover is approximately 

N co = T/{hw) In {{h 2 T/ (mg 2 )) 1 / 3 ) = T/(3hu) ln(t/2). (49) 

Since t 1 / 3 1 (see text below Eq. (l24l)). this equation can be inverted to give a crossover 
temperature 

T co = Nhw/\n ({Ntfu/img 2 )) 1 '*) . (50) 

A comparison of this formula with a numerical calculation using Yang- Yang thermodynamics 
shows very good agreement [45]. Since t 3> 1 at the crossover, Eq. ( |49l shows that the ratio 
N co /Nd, where Nd = fkj/T is the atom number at degeneracy, is larger than one at the crossover. 
Thus, even considering the extensive variable N, the degenerate ideal gas regime is in principle 
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identifiable. However, the ratio N co /Nd only grows as logarithm of t and it is in practice difficult 
to have N co /N d very large. 

3.3.2 Validity of the local density approximation 

All the previous results use the local density approximation, which requires that the correlation 
length l c of the gas be much smaller than the scale L of variation of the density. At the crossover, 
the correlation length of the gas is about l c ~ £ = %/ ^Jnign, as seen in section[2l To estimate L, 
let us approach the crossover from the ideal gas regime. The density profile of the central part of 
the cloud, obtained using Eq. (fT6b and the local chemical potential fj,(z) — fj, — muPz 2 /2, turns 
out to be a lorentzian of width ^J\^\/muj 2 . Thus, L ~ ^y\fi \/muj 2 ~ (gTco/mhw 3 ) 1 ^ 3 at the 
crossover. We thus find that the condition of validity of the local density approximation, l c <C L, 
can be rewritten as 

lu<^uj co = — 3 — = — , (51) 

a n 

a result which has been derived in [45]. 

If the local density approximation (I5TT) is not satisfied, the discrete structure of the trap energy 
levels has to be taken into account. In the opposite limit, u ^> oj co , the quantization of energy 
levels plays a role while the gas is still described by an ideal Bose gas. Then, it has been shown 
in [53] that one expects a condensation phenomenon to occur at a temperature 

T c = Nfiw/\n(2N). (52) 

In contrast to the crossover described in the previous section (referred to now as the interaction- 
induced crossover), this is a finite size phenomenon since T goes to zero when the trap confine- 
ment uj goes to 0, Nuj being fixed. This condensation phenomenon will dominate the interaction 
induced crossover when Tc > T co . This condition is equivalent to uj 3> uj co , which shows 
consistency of our analysis. 

Experimentally, the condition (|5TT) to observe the interaction induced crossover is very easily 
satisfied: using Eq. (l42l) . the condition (|5"TI) reduces to 

uj < u ± (T/hu ± ) 2 ^(a/l ± ) 2 ^. (53) 

One can check that, for most alkali atoms, in trapping potentials with u± ranging from 1 to 
several tens of kilohertz and for temperatures between 0.1Huj± and fko±, this condition is easily 
fulfilled, unless a is extremely small (a < 0.1 nm). Thus, one expects that a trapped ID gas 
undergoes the interaction induced crossover towards a quasi-condensate and that the local density 
approximation is valid to describe the gas. 

3.4 3D physics versus ID physics 

Experimentally, one expect a crossover from a one-dimensional behavior to a three dimensional 
behavior as the temperature of the gas increases and, at large enough temperature, one expects 
to recover the physics of a three-dimensional gas. The physics of a 3D gas is very different from 
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that of a ID gas. The most striking difference is that, even in the absence of interactions, a 3D 
Bose gas undergoes a phase transition towards a BEC due to saturation of the population of the 
excited states. This is in contrast to ID gases where, in absence of interactions between atoms, 
the gas behaves, for any density, as a thermal gas in which bosonic bunching is present. For 
weakly interacting gases, in both ID and 3D gases, a transition towards a (quasi-)condensate is 
expected. However, these transitions are different in nature and this difference can be captured 
by studying the validity of mean field theories in both cases. 

In 3D weakly interacting gases (pa 3 <C 1), the effect of interactions between atoms on the 
onset of Bose Einstein condenstation is very small. This is why 3D Bose gases with weak in- 
teractions are well described by mean field theories. For instance, the thermodynamics is given 
with a very good approximation by the Hartree-Fock-Bogoliubov self consistent theory [54,55]. 
In such a theory, at temperatures larger than the critical condensation temprature, the gas is 
described by the Hartree-Fock approach, in which correlations between atoms introduced by 
interactions are neglected. Condensation is then due, as for an ideal Bose gas, when the den- 
sity reaches 2.&12.../ X dB . For higher densities, a non zero condensate wave function appears, 
which is the order parameter of this second-order phase transition. The experimental value of the 
critical temperature in weakly interacting ultra-cold Bose gases is in good agreeement with this 
theory [56]. 

However, even for weakly interacting gases, such a mean field theory is expected to fail very 
close to the critical point of temperature T c . This is due to the large long wave length fluctuations 
that develop in the vicinity of the transition. In the condensate side, i.e. for T < T c , the Hartree- 
Fock-Bogoliubov self consistent theory is valid only if the fluctuations of the condensate wave 
function, averaged over a volume of the order of the correlation length, are smaller than its 
mean-field value. This is the so-called Ginzburg criteria and it gives [55] 



T, 



> ap 1/A . (54) 



c 



The same criterion (up to an absolute value) is true above T c . The region around the transition 
where \T C — T\/T c is of the order or smaller than ap 1 ^ 3 is not expected to be described by a 
mean-field. Beyond mean-field effect include a modification of the transition temperature. Since 
interactions tend to decrease long wave length density fluctuations, they favor the appearance of a 
condensate and, for small paramater ap 1 / 3 , an increase of the critical temperature is expected [57- 
59]. Such a modification is very small in cold atom experiments and has never been observed. 
A second non mean-field effect is the modification of the critical exponent that describes the 
divergence of the correlation length in the vicinity of the critical point. Measuring beatnodes 
between the atomic field extracted at different places in the atomic cloud, the critical exponent 
was measured recently in dilute atomic gases, in agreement with beyond mean-field theories [60]. 

The physics is very different in ID systems, since long wavelength fluctuations play an en- 
hanced role compared to 3D systems. The crossover towards a quasi-condensate is, in ID gases, 
a phenomenon driven by interactions. More precisely, the crossover towards a quasi-condensate 
is produced by the correlations between atoms brought by the interactions. It cannot be captured 
by the Hartree-Fock theory because Hartree-Fock theory neglects correlations between atoms in- 
troduced by interactions. Thus, in real systems the failure of the Hartree-Fock theory to describe 
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the appearance of a (quasi-)condensate is a signature of the ID nature of the physics involved. 

4 Experiments 

In this section we will discuss several experiments that have been carried out in both Orsay and 
Amsterdam using atom chips which probe the ideas discussed in the previous sections. Atom 
chip setups are very well suited to study one-dimensional geometry since very tight atom guides 
are easily realised by going close to a current-carrying micro-wire. The atom chips which were 
used in the experiment presented below are sufficiently similar that we will attempt to describe 
both at once. We will refer the reader to the individual experiments for more detailed information. 
The atom chips we used employed current carrying wires to create magnetic trapping fields for 
87 Rb atoms in the F = 2,m F = 2 state. Magneto-optical traps, laser cooling and evaporative 
cooling were used to load atoms into the chip-based traps, which tended to be highly confining 
but rather shallow. Typical currents were on the order of a few amperes and the atoms were 
at a distance of several tens of microns from the wire surface. Typical transverse confinement 
frequences (u±/2n) were about 3 kHz, while longitudinal frequencies were on the order of 10 
Hz. This transverse frequency corresponds to a temperature fauj±/kB of 144 nK, and evaporative 
cooling was able to reach a temperature equal to or slightly above this value. For Rb atoms, with 
3D scattering length a = 5.24 nm, the energy scale E g corresponds to 0.20 nK, and 144 nK in 
reduced temperature units corresponds to t = 720. Since the longitudinal trapping potential is 
roughly harmonic, the linear atom density varied in space, and thus a single sample permits one to 
probe a large range in density at constant temperature. A single density profile thus corresponded 
to a horizontal line in Fig.Q] The value of the parameter 7 was typically between 1CT 1 to 1CT 3 . 
The data consisted of absorption images of the cloud, taken either in situ or after a very short 
expansion time. Temperature measurements were made by fitting the wings of the cloud, or by 
fitting to the Yang- Yang model (see description below). 

The first set of measurements we describe are simple observations of the density profiles of 
nearly one-dimensional gases on an atom chip. The measurements were carried out with two pur- 
poses in mind. In the first measurements, carried out in Orsay, emphasis was placed on proving 
that in the region of the crossover between the ideal gas and quasi-condensate regimes, no theo- 
retical approach which neglected interaction induced correlations between particles, in particular 
the Hartree-Fock approach, could explain the profiles. In the second set, carried out in Amster- 
dam, it was shown that the exact thermodynamic treatment accounted very well for the entire 
observed profile, notably when the gas was in the crossover regime. After examining the profiles 
we move to another type of measurement in which the absorption images were analyzed to give 
information about density fluctuations. Although these measurements where chronologically the 
first, we will treat them last. 

4.1 Failure of the Hartree-Fock model 

A typical density profile is shown in Fig. @] Superimposed on the data are three different theo- 
retical predictions. The dashed line shows the profile as predicted by the ideal gas model. In the 
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wings of the profile this model should be valid, and indeed the fit to the wings of the distribution 
was used to deduce the temperature and the chemical potential of the gas. Clearly, however the 
ideal gas prediction begins to rapidly deviate from the data, because, without interactions, a ID 
Bose gas can accommodate arbitrarily high densities at a given temperature. The dash-dotted 
line shows the prediction of the quasi-condensate model Eq. (l44j) . at the same chemical potential 
as was found by fitting the wings. This model accurately reproduces the high density part of the 
distribution, but not the presence of so many atoms in the wings of the distribution. 

The Hartree-Fock theory is a variational method in which the atoms are described by a gas 
of non interacting bosons subject to an effective potential Vhf due to the mean field of the other 
atoms. Minimizing the free-energy of the gas, one finds 

V HF (v) = 2g 3D p (55) 

where g-^D = <i7ih 2 a/m is the 3D coupling constant and p is the 3D gas density. This theory 
is thus self consistent, since for a given chemical potential and temperature, p depends on Vhf- 
The factor 2 reflects the bunching, which is present in the Hartree-Fock approximation since the 
gas is descibed by a gas of non interacting bosons. 

Using minimisation techniques, the Hartree-Fock density profile was calculated in Ref. [6] 
for the experimental three-dimensional trapping potential and for the temperature and chemical 
potential found by fitting the wings of the distribution. One sees that the Hartree-Fock density 
profile, shown as a solid line in Fig. HI reproduces the wings of the density profile, and does 
not diverge as does the ideal gas profile. It does not however, reproduce the high density part 
of the profile. Moreover, the Hartree-Fock calculation shows that the Hartree-Fock gas is far 
from being saturated: the population of the ground state is very small, and no condensation is 
expected according to this mean-field model. The excess of atoms in the center is the onset of a 
quasi-condensate, although the cloud is not deep into the quasi-condensate regime. This peak is 
formed by the effects of interactions altering the two body correlation function so as to lower the 
interaction energy relative to a Hartree-Fock gas at the same density. 

4.2 Yang- Yang analysis 

Two more examples of axial density profiles measured [8] at two different temperatures and a 
peak linear density of « 50 pmr 1 are shown in Fig. [5] These profiles were fit to the model 
based on the exact Yang-Yang solutions described in Sec. 13.21 The fits are shown in the Figure as 
continuous curves, and the resulting temperature T and chemical potential fi are also indicated. 
The chemical potential p and the temperature T are the only free parameters in the model, and 
it was found that the full set of in situ measurements could be explained by the Yang- Yang- 
based model [8]. For comparison, The ideal-gas prediction and the quasi-condensate prediction 
are also shown. Clearly, the Yang- Yang-based model describes the entire profiles well, while 
the approximate models fail, in particular the smooth crossover between the two approximate 
models in the region where fx(z) ~ is captured very well by the model. 

This analysis was further corroborated by measurements of the axial momentum distribution 
[8], obtained using Bose gas focusing [37]. The tails of the momentum distribution were used to 
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Figure 4: Failure of Hartree-Fock theory in a quasi-lD gas. The experimental profile (crosses) 
is compared with the profiles expected for a quasi-condensate ( dotted-dashed), for an ideal Bose 
gas (dashed), and to the profile predicted by the Hartree-Fock theory (continuous line) for the 
same temperature and chemical potential. The vertical axis is the number of detected atoms per 
6 jum longitudinal pixel. The temperature of the gas was T = 360 nK= 2.75faj±. Adapted from 
Ref. [6]. 



extract temperatures, and these were found to agree very well with the temperatures derived from 
the Yang- Yang fit to the in situ data. The full momentum distribution is not a thermodynamic 
quantity, and can thus not be obtained directly from the Yang- Yang analysis. 

The similarity of the measured density profiles of Figures Ob) and 0] clearly suggests that 
the same physics of an interaction-induced crossover applies to both experiments. Although it is 
tempting to apply the Yang-Yang-based analysis of Sec . 13 .21 also to the data of Fig.0], this has not 
been done. It is likely that the result would not be quantitatively accurate, because at the higher 
linear densities and temperatures of Fig. HI the validity limits of the model of Sec. l3.2l are reached 
near the peak of the profile (since both fi huj± and T « fko±). Interactions are then expected 
to also play a role in the radially excited states, and also interactions among the different radial 
states will be significant. 

4.3 Measurements of density fluctuations 

As we have emphasized in Sec. [21 the transition towards a quasi-condensate in ID gases is 
characterized by the inhibition of atom bunching, the large density fluctuations characteristic of 
a thermal Bose gas. A direct measurement of the density fluctuations through the crossover thus 
captures an essential characteristic of the crossover. 

The measurement of density fluctuations proceeds similarly to the density profile measure- 
ments. The difference is that many (about 300) profiles are acquired and, roughly speaking, 
for each observation pixel, we compute the variance of the density measurements as well as the 
mean. We can relate this variance to the density fluctuations predicted by various theoretical 
approaches as described below. 
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Figure 5: Comparison of experiment to Yang-Yang thermodynamics. The model described in 
Sec, \3.2\ is fit to two examples of measured linear densities rii (dots) in the Amsterdam experiment 
[8, 38]. The resulting fits (continuous curves) yield chemical potential fi and temperature T as 
indicated. Dotted curves: ideal-gas profile at the same temperature and chemical potential 
exhibiting divergence for = 0. Dashed curve in (b): quasi-condensate profile with the 
same peak density as the experimental data. In these experiments ftw± = 158 nK. Adapted from 
Ref. [8J. 
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The measurements involve several subtleties requiring careful normalizations and corrections 
of the data. These are described in detail in [7,61]. The measurement requires a high degree of 
reproducibility in the data. The atom chip geometry permits the construction of a very compact 
apparatus with low sensitivity to vibration. The images contain not only noise due to atom 
fluctuations, but also photon shot noise. The photon noise must be carefully characterized and 
subtracted. Examples of the data are shown in Figs. [6] and [7J 

4.3.1 A local density analysis 

The pixel size in the experiment is A = 6 fim. The pixel size is much larger than the correlation 
length of the gas which is always smaller than a micron in these experiments, but much smaller 
than the longitudinal length scale of mean density variation. Thus, the data should reproduce 
number fluctuations predicted in a longitudinal local density treatment. More precisely, the gas 
contained in the pixel located at position z can be described as a gas, confined transversely by the 
transverse potential of frequency cu± and confined longitudinally by a box like potential of size 
A. The properties of this slice, which can exchange energy and particles with the rest of the gas, 
is well described within the grand-canonical ensemble. The energy shift V(z) of this slice can 
be converted to a shift — V(z) of the chemical potential. This is the local density approximation, 
already discussed in sec. [3j Since A is large compared to correlation length of the gas, the 
boundary conditions used to compute thermodynamic quantities are all equivalent and we use 
the periodic boundary conditions in the following. 

Within the local density approximation, the confinement potential V(z) is irrelevant to an- 
alyze the atom-number fluctuations. The atom number fluctuation SN 2 in each pixel depends 
only on the temperature T and on the local chemical potential. Equivalently, 5N 2 is a function 
of T and (N), since the linear density is a monotonically increasing function of the chemical 
potential. We thus choose, for each cloud temperature, to represent the measured atom number 
fluctuation as a function of the mean atom number in the pixel. Experimental results are shown 
in Fig.[6]and fig.[7J 

4.3.2 Ideal gas regime : observation of bunching 

If the gas within a pixel can be considered ideal, we can use the results of Sec. 12.21 The fluctua- 
tions of atom number rij in each one-atom quantum state |i) are 

(n 2 ) - (m) 2 = (m) + (n,) 2 . (56) 

The fluctuations of the total atom number N are thus 

(N 2 )-(N) 2 = (N) + J2(^) 2 , (5V) 

i 

where the sum is performed over all the quantum states. The mean values (n^) are given by the 
Bose distribution and the fluctuations of N are easily computed. 
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Figure 6: Density fluctuations of a gas on an atom chip. The atom number variance in an 
observation pixel (6 \im) is plotted as a function of the mean number. The open cirles are the 
fluctuations measured for a hot cloud (T — 1.3 fiK corresponds to 10Juo±) for which bunching is 
unobservable because of the large number of transverse states involved. The variance is due to 
atom shot noise. Full circles correspond to a colder cloud, at a temperature T = 2.1hw±. The 
increase in fluctuations is due to bunching. The theoretical prediction for an ideal Bose gas at 
the same temperature is given by the dashed curve. The prediction for a nondegenerate cloud, 
Eq. rf671) . is shown as the dotted curve. The degeneracy of the gas is evident. Adapted from 
Ref. [7]. 



400 
300 
200 
100 




oT~ 1.3 pK 
• T~ 190 nK 




i i i 



o 

o 







100 



200 



300 



400 



N 



Figure 7: Density fluctuations in the quasi- condensate regime. The dashed-dotted curve is the 
prediction for an ideal Bose gas at the same temperature as in figure^ The dashed curve is the 
prediction for a quasi-condensate. In units of transverse energy, the temperature is T = lAhuj±. 
Adapted from Ref. [7]. 
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A rough calculation is as follows: if M quantum states are populated with similar popula- 
tions, Eq. (1571) simplifies to 



(N 2 )-(N) 2 = (N) + (N) { -^. (58) 

The first term of the right hand side is the shot noise term, expected for uncorrelated, statistically 
independent atoms. The second term on the right hand side is the effect of the bunching. We 
see from this expression that as long as (N)/M is much smaller than 1, the bunching term is 
negligible compared to the shot noise term. The ratio (N)/M is approximately the phase space 
density of the gas and is much smaller than 1 if the gas is non degenerate. Thus one expects 
the measured atom number fluctuations to be dominated by the shot noise term for gases at high 
temperature. This is observed experimentally for non-degenerate clouds, as shown by the open 
circles in Fig. [6] The linearity of the measured value of (iV 2 ) — (iV) 2 versus (N) shows that the 
fluctuations are given by the shot noise. The fact that the slope is smaller than the expected slope 
of 1 is due to the fact that the optical resolution (about 10 /im) is larger than the pixel size [7]. 

One can also give a more precise calculation of the fluctuations. For this purpose, we in- 
dex the quantum states by the integer n x and n y , which label the transverse vibrational levels, 
and the longitudinal wave vector k z , which takes values in multiples of 2ix / A. For a highly 
nondegenerate gas, ^> T, the population of each state is given by the Boltzmann law 

fan^O = Ae-^/^^< n - +n »»/ T , (59) 
where the normalization factor A is 

A=^^(l-e^ T ) 2 . (60) 



Inserting this into Eq. (1571) . we obtain 

(iV 2 ) - (iV) 2 = (N) + (iV) 2 4^tanh 2 (^ ± /2T). (61) 

v 2 A 

We thus recover an expression similar to Eq. ( |58l ), with M = \f2A/(\ dB tanh 2 (ftc<;_|_/2T)). The 
tanh term accounts for the number of populated transverse states. The term y/2 A /XdB, which 
accounts for the longitudinal states, may be recovered by a semiclassical analysis: the volume 
of the occupied phase space is Vt ~ AymT and the number of quantum states contained in this 
volume is of the order of VL/h. Equation (I6TI) is valid if the gas is non degenerate. When the gas 
becomes degenerate, the distribution of the mean occupation number (n) versus the state energy 
becomes more peaked around zero. This amounts to a reduction of the effective number of 
occupied states M and the effect of bunching is larger than the prediction of Eq. (|6Tb . For highly 
degenerate gases, Eq. (I6TT) underestimates the true fluctuations, which become large compared 
to the shot noise level. 

The bunching effect is quite clear for a cold enough cloud as shown in Fig. [6j In this ex- 
periment, the bunching term is even larger than the shot noise term, indicating that the gas is 
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degenerate. The degeneracy is also shown by a comparison of the data with Eq. (|6TT) shown 
as a dotted line. This equation, valid for a non degenerate gas, underestimates the measured 
fluctuations. On the other hand, a calculation of Eq. (1571) using the true Bose occupation factor 
is in much better agreement with the data. This comparison shows that, at least as concerns 
fluctuations, the gas is well described by an ideal, degenerate Bose gas. 

4.3.3 Quasi-condensate regime: saturation of atom number fluctuations 

At sufficiently high density and low temperature, repulsive interactions between atoms are no 
longer negligible. As described in section [2} one expects the interactions to reduce the density 
fluctuations to lower the interaction energy. The gas then enters the quasi-condensate regime. 
For the temperature T = 2.1 fko± of the data in Fig.[6l using Eq. (|42j) and assuming a purely ID 
gas, Eq. d24l) gives a density at the crossover of about 130 atoms per pixel. Although Eq. (l24l) 
does not apply since the gas is not purely ID, this rough estimate shows that the crossover 
to a quasi-condensate is achievable at slightly higher atom number and/or lower temperature. 
Measurements of atom number fluctuations in a regime where the cloud center is in the quasi- 
condensate regime are shown in Fig. [7] Whereas at low atomic density the measured fluctuations 
are in agreement with the ideal Bose gas prediction, one sees a saturation of the fluctuations at 
higher densities. 

To calculate the fluctuations, we first suppose the gas to be purely one dimensional, with 
a coupling constant g given by Eq. (|42"1) . In a local density approximation, we consider atom 
number fluctuations in a longitudinal box of length A in equilibrium with a reservoir of energy 
at temperature T and a reservoir of particles at chemical potential fi. As explained in sec. l2.3[ the 
Hamiltonian is quadratic in Sn in the quasi-condensate approximation and Sn can be expanded 
as a sum of independent modes indexed by the wave vector k. The atom number fluctuations 
N — (N) are obtained by integrating Sn over the pixel size. Thus, the only excitation that leads 
to atom number fluctuation is the zero momentum mode. Its energy, derived from Eq. (1321) , is 

H k=0 = A 9 -5nl (62) 

Using the equipartition theorem, we find Sri 2 , = T/(gA). The atom number fluctuations, which 
are (ASn ) 2 , are thus given by: 

AT 

(N 2 ) - (N) 2 = — . (63) 
9 

Thus, we expect the atom number fluctuations to be independent of (N). The shot noise term 
is not present in this quasi-condensate regime: interactions between atoms prevent even the 
shot noise fluctuations and at temperature smaller than gn, one expects to observe sub-shotnoise 
fluctuations. This feature is also seen in Fig. |U where ^ 2 )(0) goes below unity in the exact 
solution. 

In the experimental results shown in Fig. [71 the typical interaction energy is about 0.7huj±. In 
these conditions, the transverse degrees of freedom cannot be neglected and the result of Eq. (1631) 
must be corrected. More precisely, the phonons, which are longitudinal density waves, are as- 
sociated with a breathing of the transverse shape of the cloud. For phonons of frequency much 
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smaller than the transverse frequency, the transverse shape of the cloud follows adiabatically the 
ground state equilibrium state for the local linear density. We denote by E eq {n) the energy of 
the gas per unit length for a linear density n. The phonon Hamiltonian of Eq. (l32l) is the term 
of the Hamiltonian of order two in 5n and V6. Thus, the interaction term of the phonons is 
E- mt = \L(d 2 E eq / dn 2 )8n\. Since dE eq /dn is the chemical potential of the gas (to zero order 
in 5n), we can rewrite the former expression as E int = ^L(dfj,/ 'dn)8n\. In particular, the zero 
momentum term is 

Hk=o = \^nl (64) 
Then, the atom number fluctuations are 

Although we derived this expression in the approximate quasi-condensate theory using expansion 
of the Hamiltonian to second order in 8n, we recover here a well known result of statistical 
physics. More precisely, as shown in [62], Eq. (|65l) holds for any system in equilibrium with a 
particle reservoir at chemical potential ji and with an energy reservoir of temperature T. 

To apply Eq. (1651) to the experiment, we need the equation of state /u(n). Using Eq. (1441) and 
Eq. (|65l ), one can compute the expected fluctuations. Fig. [7] shows that the results are in fairly 
good agreement with the measured atom number fluctuations. 



5 Conclusion 

We hope that we have given the reader a useful overview of the physics of ID gases in the weakly 
interacting regime. These systems are rich and manifest several different regimes separated by 
smooth crossovers. It is often necessary to appeal to many different physical models to under- 
stand them. The existence of exact solutions allows us to test the models and to explore their 
validity in the face of highly non-trivial many body correlations. Although the experimental and 
theoretical work we have described is quite extensive, we believe that much work remains to be 
done. The study of fluctuation phenomena is still at an early stage. For example the exact ther- 
modynamics should admit a careful comparison with data such as that in Fig. [Vj and improved 
experiments should probe larger parameter ranges and possibly even permit measurements of 
the correlation length. It may also be possible in the near future to enter the strongly interacting 
regime using an atom chip. Experiments are also capable of measuring momentum distribu- 
tions [8], but so far no quantitative theoretical comparison has been made. Finally measurements 
of fluctuations and correlations in momentum space are also experimentally feasible [32]. 
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